Illustrating the Geometry of Coherently Controlled 
Unital Open Quantum Systems 

Corey O'Meara, Gunther Dirr, and Thomas Schulte-Herbriiggen* 

dated: Aug. 15. 2011 



Abstract — We extend standard Markovian open quantum sys- 
tems (quantum channels) by allowing for Hamiltonian controls 
and elucidate their geometry in terms of Lie semigroups. For 
standard dissipative interactions with the environment and dif- 
ferent coherent controls, we particularly specify the tangent 
cones (Lie wedges) of the respective Lie semigroups of quantum 
channels. These cones are the counterpart of the infinitesimal 
generator of a single one-parameter semigroup. They comprise 
all directions the underlying open quantum system can be 
steered to and thus give insight into the geometry of controlled 
open quantum dynamics. Such a differential characterisation is 
highly valuable for approximating reachable sets of given initial 
quantum states in a plethora of experimental implementations. 
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I. Introduction 

Extending quantum channels by allowing for Hamiltonian 
control turns them into interesting and important examples of 
geometric control of open systems. While for closed systems, 
the theory of Lie groups provides a rich structure to address 
questions of reachability, accessibility, and controllability [1], 
already simple open quantum systems come with the intricate 
geometry of Lie semigroups [2,3]. For instance, in most 
closed systems the reachable set to an initial state po simply 
is the orbit O G (p ) := {Gp a G~ 1 \G e G} of a unitary 
subgroup G whose Lie algebra can be identified easily via Lie 
closure, while in open systems reachable sets are much more 
difficult to determine explicitly. Thus in view of controlling 
open quantum dynamics, in [4] we systematically related the 
framework of completely positive semigroups [5-11], which 
is well established in quantum physics, with the more recent 
mathematical theory of Lie semigroups. An early example 
confined to single-qubit systems can be found in [12]. 

More precisely, for exploiting the power of systems and con- 
trol theory in open quantum dynamics, the system parameters 
have to be characterised first, e.g., by input-output relations 
in the sense of quantum process tomography. The decision 
problem whether the dynamics of the quantum system thus 
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specified is Markovian to good approximation has recently 
been analysed [13, 14]. Moreover (time-dependent) Markovian 
quantum channels were elucidated from the viewpoint of 
divisibility [13] thus paving the way to Lie semigroups [4]. 
Following up, this work sets out to determine the geometry of 
quantum channel semigroups in terms of their tangent cones 
(Lie wedges) for a number of coherently controlled standard 
unital channels in a unified frame in line with [15]. 

For the first time, here we explicitly parameterize the set of 
all possible directions an open quantum system under coherent 
controls may take — its Lie wedge. Thereby, we heavily 
exploit the fact that the set of all reachable quantum maps 
governed by a controlled Markovian master equation consti- 
tutes a Lie semigroup [4]. Previous characterizations of reach- 
able sets for unital open quantum systems by majorization 
techniques, e.g., [16], become increasingly inaccurate once 
full controllability of the Hamiltonian part (condition (H) vide 
infra) is violated, which for growing number of qubits happens 
in all experimentally realistic settings. In contrast, the Lie- 
semigroup tools presented here do not require condition (H) 
and carry over to multi-qubit systems without the draw-back 
of increasing inaccuracy. 

II. Theory and Background 

We start out by recalling some basic notions and notations 
of Lie subsemigroups [2] and their application for character- 
ising reachable sets of quantum control systems modelled by 
Lindblad-Kossakowski master equations [4]. 

A. Lie Semigroups 

To begin with, let G be a matrix Lie group, i.e. a 
group which is (isomorphic to) a path-connected subgroup 
of GL(n,R) or GL(n,C) for some n G N, and let g be 
its corresponding matrix Lie algebra. Thus g is (isomorphic 
to) a Lie subalgebra of g[(n,R) or g[(n,C). Then a subset 
S C G which is closed under the group operation in the sense 
S • S C S and which contains the identity 1 is said to be a 
subsemigroup of G. The largest subgroup within S is written 
E(S) :=SnS _1 . 

Furthermore, a closed convex cone ro C g is called a wedge. 
The largest linear subspace of ro is denoted E(n>) := ron(— ro) 
and it is termed the edge of the wedge ro. Now, ro G g is a 
Lie wedge of g if it is invariant under the adjoint action of the 
subgroup generated by the edge E(to), i.e. if it satisfies 

e A roe" A = ro (1) 

(or equivalently e adA (ro) = ro) for all A G E(n>). Note that 
the edge of a Lie wedge always forms a Lie subalgebra of g. 

Moreover, for any closed subsemigroup S of G we define 
its tangent cone L(S) at the identity 1 by 

L(S) := {A G g | cxp(tA) G S for all t > 0} . (2) 

Then one can show that L(S) is a Lie wedge of g satisfying 
the identity E(L(S)) = L(E(S)). Yet, the 'local-to-global' 
correspondence between Lie wedges and closed connected 
subsemigroups is much more subtle than the correspondence 
between Lie (sub)algebras and Lie (sub)groups: for instance, 



several connected subsemigroups may share the same Lie 
wedge ro in the sense that L(S) = L(S') for S ^ S', 
or conversely there may be Lie wedges ro which do not 
correspond to any subsemigroup, i.e. ro = L(S) fails for all 
subsemigroups S G G. 

Therefore, one introduces the important notion of a Lie 
subsemigroup S which is characterised by the equality 

S = (expL(S)) s , (3) 

where the closure is taken in G and (expi(S))s denotes 
the subsemigroup generated by expL(S), i.e. (expi(S))s := 
{e Al ---e An \n G N, A t , . . . , A n G L(S)}. Moreover, a 
Lie wedge ro is said to be global in G, if there is a Lie 
subsemigroup S C G such that 

L(S) = fo • (4) 

Thus, one has the identity S = (cxpro) s . 

Whenever a Lie wedge ro C g specialises to be compatible 
with the Baker-Campbell-Hausdorff (BCH) multiplication 

A * B := A + B + | [A, B] + ■ ■ ■ = log(e A e B ) VA, Bern 

(5) 

defined via the BCH series, it is termed Lie semialgebra. For 
this to be the case, there has to be an open BCH neighbourhood 
B G Q of the origin in g such that (ro n B) * (ro H B) G ro. 
An equivalent definition for being a Lie semialgebra is given 
by the tangential condition 

[A, Ta tn] c T A xo for all A G to , (6) 

where T^ro denotes the tangent space of ro at A defined by 

TaXo := (A 1 - n to*) ± . (7) 

Here A 1 - denotes the orthogonal complement of A and 
ro* := {A G g | (A, B) > for all B G to} the dual 
wedge — both taken with respect to the standard trace inner 
product. The conceptual importance of Lie semialgebras roots 
in the fact that — in Lie semialgebras — the exponential map 
of a zero-neighbourhood in L(S) yields a 1-neighbourhood 
in S. In contrast, as soon as ro is merely a Lie wedge that 
fails to carry the stronger structure of a Lie semialgebra, 
there will be elements in S that are arbitrary close to the 
identity without belonging to any one-parameter semigroup 
completely contained in S. For more details and a variety 
of illustrative examples, we recommend [2] and [17], where 
the respective introduction does provide a lucid overview of 
the entire subject. The connection between Lie semialgebras 
and time-independent Markovian quantum channels has been 
worked out in detail in [4]. 

With these stipulations, the frame is set to describe the 
time evolution of Markovian (i.e., memory-less) open quantum 
systems in the differential geometric picture of Lie wedges. 

B. Markovian Quantum Dynamics and Quantum Channels 

Markovian quantum dynamics is conveniently described by 
a linear autonomous differential equation 

X(t) = -£X(t), (8) 
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where X(t) usually denotes the state of a quantum system 
represented by its density operator p(t), i.e. p(t) = p(ty , 
pit) > 0, and tr p(t) = 1. Here and henceforth, (•)' denotes 
the adjoint (complex-conjugate transpose). For ensuring com- 
plete positivity, C has to be of Lindblad form [10], i.e. 

C(p) = iad H (p)+T L (p) , (9) 

with ads, (p) '■= [Hj,p] and 

Fl(p) := | V k V *P + P V kVk ~ 2V k pV2 , (10) 

k 

Here, the Hamiltonian H is assumed to be a Hermitian N x N 
matrix while the Lindblad generators {Vk} may be arbitrary 
N x N matrices. The resulting equation of motion (8) acts on 
the vector space of all Hermitian operators, her(iV), and more 
precisely, leaves the set of all density operators posi(N) := 
{p E t)tt(N) \ p — p^,p> 0,trp = l} invariant. 

In [4] it was shown that the set of all Lindblad generators 
{— £} has an interpretation as a particular Lie wedge. To see 
this, consider the group lift of (8), i.e. now X(t) denotes an 
element in the general linear group GL(f)ez(N)). Moreover, 
define the set of all completely positive (cp), trace-preserving 
invertible linear operators acting on het(iV) as P cp , i.e. 

pep ._ jji g GL(t)cx(N)) | T is cp and trace-preserving} 

and let Pq P denote its connected component of the identity. 
Then, P cp is exactly the set of so-called invertible quantum 
channels. A quantum channel T is said to be time independent 
Markovian or briefly Markovian, if it is a solution of (8). Thus 
T = c- tc for some fixed Lindblad generators C and some 
t > 0. Furthermore, T is time dependent Markovian if it is 
a solution of (8), where now C = C(t) may vary in time 
(for terminology see also [13,14]). Finally, we will denote 
the set of all time independent Markovian and time depen- 
dent Markovian quantum channels by MQC and TMQC, 
respectively. Then, with regard to the work by Lindblad [10] 
and Kossakowski [9], one obtains the following result [4]: 

(a) The global Lie wedge of Pq P is given by the set of all 
Lindblad generators of the form 

-C:=-(i&d H +T L ) (11) 

with H e het(iV) and T L as in (10). 

(b) The Lie semigroup 

NWI S (i2) 

clearly contains MQC and moreover it exactly coin- 
cides with the closure of TMQC thus excluding the 
non-Markovian ones in Pq P , which is most remarkable. 
While assertion (a) reformulates previous results by Lindblad 
and Kossakowski [6,9,10], part (b) is noteworthy as it also 
says that Pq P is not a Lie subsemigroup of GL(\)tt(N)). 

C. Coherently controlled Master Equations 

Controlled Markovian quantum dynamics is appropriately 
addressed as right-invariant bilinear control system [4, 18-20] 

p(t) = -C u[t) {p(t)) , P (0) epo Sl (N) , (13) 



where C u now depends on some control variable ii£l m . 

Here, we focus on coherently controlled open systems. This 
means that C u has the following special from 

C u {p) = -i&d Hu {p) -r L (p) with (14) 

m 

adf/ u := &d Hd +^uj &d H] ■ (15) 

3=1 

Note that the control terms i ad# with control Hamiltonians 
Hj £ her(iV) are usually switched by piecewise constant 
control amplitudes Uj(t) G M. The drift term of (14) is 
then composed of two parts, (i) the term i &dn d (in abuse 
of language sometimes called 'Hamiltonian' drift) accounting 
for the coherent time evolution and (ii) a dissipative Lindblad 
part r^. So C u denotes the coherently controlled Lindbladian. 
As in the uncontrolled case, system (13) acts on the vector 
space of all Hermitian operators leaving the set of all density 
operators invariant. Equivalently, one can regard (13) as an 
affine system on hcr (iV) := {H e ljcz(N) \ txH = 0}. 

In the following, we further impose unitality, i.e. we assume 
Tl(T) = 0. This ensures that (13) actually yields a bilinear 
control system on het (7V) instead of an affine one. Therefore, 
it allows a group lift to GL(het (-/V)) which henceforth is 
referred to as (£), i.e. 

(E) X(t) = -C u(t) X(t) , X(0) e GL(het (AO) • (16) 

The corresponding group lift in the affine case is more involved 
[4, 19]. Now, the system semigroup Ps associated to (£) reads 

Ps = (T u (t) = cxp{-t£ u )\t> 0,u £ R m ) s (17) 

and lends itself to exemplify the notion of a Lie wedge. To 
distinguish between different notions of controllability in open 
systems, we define three algebras: the control algebra t c , the 
extended algebra id, and the system algebra s as follows 

l c := {ia,d Hj \j = l,...,m) L ie , 

t d := (i nd Hd , i ad Hj \ j = 1, . . . , m)u e , 

(18) 

s := (C u \ Uj e K}i_ie 

= (i a.d Hd +T L ,i ad Hj | j = 1, • • • , m) Ue . 

Note that s is different from td, because it contains the entire 
drift term (iad# d +Tl) for the Lie closure, while id only 
takes its Hamiltonian component iadjf d . Then (S) is said to 
fulfill condition (H), (WH), and (A), respectively, if 

(H) i c = ad su(JV) (19) 

(WH) i d = ad su (jv) while t c ^ ad 5uW (20) 

(A) s = 0[(het o (AO) (2D 

While condition (A) respects a standard construction of non- 
linear control theory [1,21] to express accessibility, conditions 
(H) and (WH) serve to characterize different types of con- 
trollability of the Hamiltonian part of (S) in the absence of 
relaxation: Condition (H) says that the Hamiltonian part is 
fully controllable even without resorting to the drift Hamilto- 
nian, whereas condition (WH) yields full controllabilty of the 
Hamiltonian part with the drift Hamiltonian being necessary. 
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We refer to the first scenario as (fully) H -controllable and to 
the second as satisfying the (WH)-condition. Generically, open 
systems (E) given by (16) meet the accessibility condition (A) 
[22,23]. 

Finally, note that via c ladH (p) = c lH pc~ lH the Lie algebra 

ad su (jv) generates the Lie group Adgj/nv) = PSU(N) here 
acting on t)cio(N) by conjugation. 

D. Computing Lie Wedges for Controlled Master Equations 

Here, the goal is to determine the (global) Lie wedge of a 
coherently controlled unital open system (E) given in terms 
of its Markovian master equation (16) of GKS-Lindblad form. 
In view of the examples worked out in detail in Sec. Ill, here 
we sketch how to approximate a Lie wedge of a controlled 
Markovian systems in two ways, (i) by an inner approximation 
and (ii) by an outer approximation thus following [3,4]. 
Moreover for unital systems, we present two results which 
guarantee that the inner approximation is global and thus 
coincides with the Lie wedge L{Py) sought for. 

Let (E) be a unital open control system as in (16) where, 
for simplicity, the system algebra s fulfills the accessibility 
condition (A). Moreover, let 

fi s := {C u \u £R m } C 0l(het o (AT)) (22) 

be the set of all directions specified by (16). The reachable 
set Reach (Sis, 1) of (E) is defined as the set of all states 
X(T), T > that can be reached from the unity X(0) = t 
under the dynamics of (E), while the controls u(t) G M. m are 
assumed to be piecewise constant functions. In general, one 
could allow for larger classes of admissible controls, such as 
locally bounded or locally integrable ones. Yet, the closure of 
the corresponding reachable sets will not differ [20,21,24]. 

Clearly, Reach (Os , 1) takes the form of a subsemigroup 
within the embedding Lie group GL(t)cv Q (N)) in the sense of 
Sec. II-A. For instance, restricting the control amplitudes {uj} 
to be piecewise constant yields the equality Reach (f2s, 1) = 
Pe. More generally, the following result holds. 

Theorem 1 ([3]): Let Ps be defined as in (17). Then 

P E = Reach 1) = Reach (L(P E ), 1) . (23) 

In particular, Ps is a Lie subsemigroup. Furthermore, £(Ps) 
is the smallest global Lie wedge containing as well as the 
largest subset S7' of gl(t)zz (N)) which satisfies the equality 

Reach = Reach . (24) 

Due to the last property, the Lie wedge L(Ps) is also called 
the Lie saturate of f2s> cf. [3,25,26]. 

Unfortunately, for an arbitrary system (E), currently no 
procedure is known to explicitly determine its global Lie 
wedge. Yet there is a straightforward strategy to compute an 
inner approximation [3,4]. It consists of the following steps: 

(1) form the smallest closed convex cone ro containing fi^; 

(2) compute the edge E(m) of the wedge and the smallest 
Lie algebra c containing ^(fo), i.e. e := (£?(ro)}Li e ; 

(3) make the wedge invariant under the Ad-action of e by 
forming the set \J Aec Ad cxpj 4(ro); 



(4) update by taking the convex hull conv {S} of the set S 
obtained in step (3); 

(5) repeat steps (2) through (4) until nothing new is added: 
the resulting final wedge roo is henceforth referred to as 
inner approximation to the global Lie wedge L(Ps). 

Now, the crucial question arises whether the inner ap- 
proximation roo is global or not. If it is global, Theorem 1 
guarantees that ro is equal to -L(Ps). Next we present two 
results which proved quite helpful to decide the globality 
problem: The first one yields a global outer approximation ro° 
of L(Ps). Combining inner and outer approximation, the Lie 
wedge L(Ps) sought for can be determined via the inclusions 

Wo C L(Pe) C ro°. (25) 

Clearly, if the outer and inner approximations coincide, one 
is done. The second one based on the so-called Principal 
Theorem of Globality from [2] (see also Appendix A) provides 
a 'direct' method for proving globality. It will be the key tool 
to show that the inner approximations given in the worked 
examples of Sees. Ill and IV are in fact global Lie wedges. 

Theorem 2 ([4]): Let (S) be a unital controlled open sys- 
tem as in (16). If there exists a pointed cone c in the set of 
all positive semidefinite operators that act on her (A^) so that 

(1) r L G c 

(2) [c, c] C ad su(Ar) 

(3) [c, ad su(JV )] C (c - c) 

(4) Ad v cAd w C c for all U G SU(N), 

then the subsemigroup associated to (E) follows the inclusion 
Ps Q ^suiN) ' CX P(~ c ) an d hence its Lie wedge obeys 
the relation L(Ps) C ad 6u ( A r) ffi(-c), i.e. ad su(A r) ©(— c) is 
a global outer approximation to L(Ps). 

Corollary 2.1: ([4,12]) Let (E) be a unital single-qubit 
system satisfying condition (H) with a generic 1 Lindblad term 
r^. Then Ps = Adsu(2) ' exp(— c), where the cone 

c := IR + conv {Ad v T L Ad uf \ U G SU(2)} 

is contained in the set of all positive semidefinite elements in 
g[(her ). Furthermore i(P s ) = ad 6u(2 ) ©(— c). 

Theorem 3: Let (E) be a unital controlled open system 
given by (16). In addition assume that (E) meets the ac- 
cessibility condition (A) and that the Lie subgroup K which 
corresponds to the control algebra i c is closed within SU(N). 
Then, to := t c © (— c) is a global Lie wedge in g[(het (2V)), 
where c := K,|conv{ Ad v (a,d iHd +T L ) Ad^ U G K} . 
Moreover, ro is the global Lie wedge of (E), i.e. ro = L(Ps). 

Proof: (Sketch) The full proof will be given elsewhere in 
a more general context. For applying the 'Principal Theorem 
of Globality' [2] (see Appendix A), the following steps have 
to be established: 

(1) The edge of ro coincides with { c . 

(2) ro is a Lie wedge in g := gI(f)et (iV)). 

'in [4] Cor. 2.1 is stated under the above genericity assumption; yet one 
can drop this additional condition. 
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(3) 



There exists a function tp : GL(t)n (N)) — > R such 
that its differential satisfies d(p(X) AX > for all X G 
GL(her (iV)) and all A e ro. 

The differential of v? fulfills dy>(l) A > for all A £ 

that step (3) is the essential one, and an appro- 
candidate for <p is given by X n> — (X, X) := 
i(B k )X{B k )) , where B 1 ,..., B N 2_ 1 is any 
orthonormal basis of t]zt (N). ■ 



(4) 



Note 
priate 



,jV -1 



As a useful tool, we add the following Corollary, which is put 
into a broader context in Appendix A: 

Corollary 2.2 ([2]): Let G be a Lie group with Lie algebra 
g and let roo Q ro be two Lie wedges in g. Provided one has 
ro \ —Wo Q It) \ —It) [or equivalently i?(ro ) = E(to) n ro ], 
then too is global in G if the following conditions are satisfied: 
(i) to is global in G; (ii) the edge of too is the Lie algebra of 
a closed Lie subgroup of G. 

Guideline through Applications 

For illustrating the power of the Lie-semigroup formalism 
by applications, we follow a two-fold route: Sec. Ill addresses 
three paradigmatic types of bilinear control systems on R 3 , 
where the control parts of the dynamics generate easy-to- 
visualise rotations in 5*0(3). Thus Sec. Ill is meant to be 
readable without any background in quantum mechanics, yet 
it directly corresponds to single-qubit systems undergoing 
relaxation as the presented examples coincide with the so- 
called coherence-vector representation of such systems [27]. 
Therefore, the results obtained in Sec. Ill can readily be 
transferred to Sec. IV, where we address quantum channels 
in the customary explicit su(2)-representation of qubits. By 
the isomorphism so(3) = su(2), the geometry in Sec. Ill thus 
illustrates key results in Sec. IV for qubit channels. 

III. Geometry of Open Systems in R 3 

In this section, we discuss three simple introductory ex- 
amples of 'open' systems, the geometry of which can be 
envisaged as rotations in M 3 concomitant to relaxation. To 
fix notations, define the following 



o o 


1 



1 

ooo 

-10 





1 





1 





(26) 



as generators of the rotations R„{0) := e " reading 



Rx(B) ■-- 



10 o 

cos(0) - sin(e) 
sin(fl) cos(6>) 



Ry(9) 



cos(S) sin(e) 

1 
- sin(0) cos(e) 



COS(V) 

sin(S) 



- sin(e) 
cos(0) 



(27) 



So we have (H x , H y , H z )\_ le = so(3) and thereby a basis for 
the skew-symmetric matrices forming the t-part in the Cartan 
decomposition £jt(3, R) = so(3) ffisrjm(3), where the p-part is 
spanned by the symmetric matrices 



Px 



"0 0" 




"0 1" 




"010" 


1 


. Py 





. Pz ■ = 


1 


10 


10 






(28) 



Table I 
Commutation Table 





fill 


E22 


E33 


Px 


Py 


Pz 


H x 





Px 


~Px 


-2A 2:! 


-Pz 


Py 


J J ii 


-Py 





Py 


Pz 


-2A31 


~Px 


H z 




-Px 





-Py 


Px 


-2A12 


define Aj, 


'■= En 


~ E JJ 











and the diagonal 3 x 3-matrices En := e^e^ for i = 1,2,3. 
Recall that the skew-symmetric t-part and a symmetric p-part 
obeying the usual commutator relations [t, t] C t, [t, p] C p 
and [p,p] C {. For later convenience, we note commutation 
relations for the above basis in Tab. I. 

A. Example 1: Corresponds to a Qubit System with Condition 
(H) Satisfied and General Relaxation Operator 

Using definitions from above, consider the control system 
in GL(3, M) given by the equation 



X = -{A + B U )X , 



(29) 



u x H x 



u y Hy shall have 



where the control term B u 
independent controls u x ,u y S R, and the drift term A := 
Hz +To is composed of a 'Hamiltonian' component, H z , and 
a relaxation component given by the matrix 



To := diag (a, b, c) 



(30) 



with relaxation-rate constants a,b,c > 0. Since (H x , H y )u e = 
so(3), system (29) satisfies in fact condition (H) in the sense 
of Sec. II-C (i.e. without resorting to the drift component H z ). 

For explicitly computing the Lie wedge of (29) we proceed 
as in Sec. II-D. For the following calculations observe that 
H x , H y , H z belong to the t-part, while To is contained in the 
p-part of 'the' Cartan decomposition of g[(3,R). 

Step (1) of the algorithm gives the initial wedge approximation 

toi :=RH x ®RH y ®(-Ci) , (31) 

(H z + Tq). In step (2) one then readily finds 



where ci := 



Eitox) 



■ H T 



' Hi, 



(32) 



so e = (H x , H y )uc = so(3). Hence, step (3) and (4) give 



ro =so(3) © R convO S o(3)(r 



(33) 



where S O(3)(To) := {er 8 T | 6 £ 50(3)} denotes the 
orthogonal orbit of To- Here we used the trivial fact that so(3) 
lso(3)" 



is Adc;nr?irinvariant. By a well-known result of Uhlmann 2 the 



convex hull of the isospectral set Oso(3)(^o) simplifies to 

coiivC9 5O(3) (r ) ={Se snm(3) | S -< T } =: A^(r ) . 

(34) 

Defining the pointed convex cone Co := Rj Ai(To), we obtain 
ttJ =so(3) © (-c ) , (35) 

2 The result mentioned is originally stated for density matrices and SU(N). 
However, the proof in [28] immediately carries over to symmetric matrices 
and SO(N) [28-30]. 
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as final inner approximation to the global Lie wedge of (29). 

Lemma 3.1: The set tt>o is a Lie wedge of fll(3, K). Its edge 
-E(roo) is given by so(3). 

Proof: It suffices to show that the edge of too is given 
by so(3). Then the invariance of too under the Ad-action 
of E(k>o) is obviously guaranteed by construction. Clearly, 
one has the inclusion so(3) C E(Wq). Conversely, let W G 
E(n> ). Then, W = A + B with A G so(3) and B G c . 
Since -W G £(«o), there exits A' e so(3) and B' e Co 
such that -W = A' + B'. Hence A + B = —(A' + B') and 
thus A + A' = -(B + B'). Since c G st)m(3), it is trivial 
that so(3) n c = {0}. Therefore, A + A' = —(B + B') = 
and hence A 1 = -A and B' = -B. Now, B, -B G c = 
Rq M(Tq). But M(T ) is contained in the set of all positive 
semidefinite matrices and therefore we conclude B = 0. Thus 
we obtain E(tv ) = so(3). ■ 

Proposition 3.1: The set too = so(3) © (— Co) is the global 
Lie wedge to control system (29). 

Proof: This is an immediate consequence of Theorem 3 
and the fact that (29) comes from a unital GKS-Lindblad 
master equation in the coherence-vector representation. ■ 

Remark 1: Alternatively to the above proof, one could 
apply Corollary 2.1, because too is a matrix representation 
of the global Lie wedge described therein. 

For general Tq = diag (a, b, c), the Lie wedge too in Exam- 
ple 1 does not carry the special structure of a Lie semialgebra, 
cf. Sec. II-A. This can be shown by choosing a suitable 
A' G too which violates the inclusion [A', X^/too] C X^/too- 
— In contrast, for the case 3 Fq = A • 1, indeed we obtain a Lie 
semialgebra, because the BCH-product A*B obviously stays 
inside too, whenever A,BG trio- Further details and proofs 
are given in Appendix B. 

B. Example 2: Corresponds to a Qubit System with Condition 
(WH) Satisfied and Control Invariant Relaxation Operator 

Consider the control system in GL(3,M) given by 



X = -(A + B U )X , 



(36) 



where the control term B u is of the form B u := uH y with 
u 6 K and the drift term A := Tq + H z is composed of a 
'Hamiltonian' part, H z , and a relaxation part 



r := 7 diag (1,0,1) 



(37) 



with 7 > 0. Since (H y ,H z )ue = so(3) the system (36) fulfills 
condition (WH) in the sense of Sec. II-C but obviously not 
condition (H). 

Now in step (1) of the inner approximation procedure to the 
global Lie wedge of (36) one finds 



roi = WLH y © (-ex) 



(38) 



3 Note that this case exactly corresponds to the (isotropic) depolarising 
quantum channel discussed in Sec. IV. 




(b) 




Figure 1 . The Lie wedge for the system of Example 2 satisfying the (WH)- 
condition is four-dimensional: it contains (a) the convex cone Co (see (42) with 
A = 1) generated by e SH y(To + H z )e~ SH y and shown in the projection 
into the subspace spanned by {H x , H z , To}; and it comprises (b) the prism- 
shaped wedge projected into the subspace spanned by {H y , H z , To}. Note 
that the edge of the wedge is spanned by the control Hamiltonian H y . Both 
parts of the figure scale with A — > oo . 



with Ci := S.q(H, + T ) whose edge is given by 
£(roi) =RH y . 



(39) 



In step (3) we include elements obtained by conjugations 
generated by edge elements identified in step (2), i.e. elements 
of the form 



c eH y ci e~ 9Hv 



X (cos(9)H z + sm(9)H x + T ) 



(40) 



for 9 G R and A > 0. By orthogonality (r |i/ 1 ,)=0 
for v = x,y, z one readily gets a Hilbert space "H := 
span {H x , H z , Tq}, in which the edge-invariant cone elements 
of (40) can be expanded using the following short-hand 



A( sin(6)H x + cos{0)H z + T ) =: A 
Then its convex hull gives the final cone 

c := Rrt cc 



" sin(0) " 




' H x ' 


cos(S) 




H, 


1 




r 





" sin(0) " 








{ 


cos(S) 




Hz 


| 6 G K J 




1 




r 





(41) 



(42) 



— a classical 3-dimensional 'ice cone', cf. Fig. 1(a). By 
construction, Co remains Ad cxp £;( 1T ,)-invariant. Finally, since 
H y _L Co, the Lie wedge itself admits the orthogonal decom- 
position 



ro 



IHy © (-C ) 



(43) 



Proposition 3.2: The set tt>o = M.H y © (— Co) is the global 
Lie wedge to control system (36). 

Proof: The Lie wedge property of roo can be derived as 
in Example 1. Then the globality of too follows again from 
Theorem 3. ■ 
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Remark 2: For clarity, let us denote the Lie wedge of 
Example 1 for Tq as in (37) by to while too still refers to 
the Lie wedge of Example 2. Clearly, one has too C to and 
E(tVo) = i?(ro )nroo. Hence globality of ro also follows by 
Corollary 2.2 and the globality of ro . 

Note that the Lie wedge too in Example 2 does not specialise 
to the form of a Lie semialgebra as can readily be verified by 
a counter example: According to (42), choose B 6 Mo as 
B = T + H x (recalling T := E n + E 33 and A = T + H z ). 
Then by the commutator relations of Tab. I, the BCH product 



A*B = 2r + H x + H, + h(H y +p x +p x ) + 



(44) 



immediately leads outside the Lie wedge too, e.g., by the non- 
vanishing component p x +p z - (NB: This argument can be made 
rigorous by introducing a scaling factor t to give tA ★ tB). 

Finally, the edge of too in Example 2 is .E(too) = M.H y (for 
7 > 0, see Fig. 1) while in the limit of a closed system, i.e. 
for 7 = 0, it turns into the entire Lie algebra so (3). 

C. Example 3: Corresponds to a Qubit System with Condition 
(WH) Satisfied and General Diagonal Relaxation Operator 

Consider the contol system in GL(3, M) given by 



X = -(A + B U )X , 
where A := T + H z , B := uH y , and 

To :=7diag(l,l,2) 



(45) 



(46) 



with ael and 7 > 0. So for approximating the correspond- 
ing Lie wedge, we take the first step to be 



roi := 



(-Ci) , 



(47) 



with Ci := Rq (H z + To) and edge given by the span of H y 
— the control 'Hamiltonian'. Again, in step (2) we identify 
the conjugation to be brought about by Adexp^n,) acting on 
the drift terms such as to give in step (3) the set 

c 3 = R+{R y (9){T + H z )R y (8) T \ 9 e M} (48) 

with the ^-component brought about by the conjugated drift 

R y (8)H z R y (0) T = cos(6)H z + am(9)H x (49) 

and the p -component reading 



R y {9)ToR y (9) 
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l+sin 2 (0) sin(0) cos(S) 
1 

sin(0)cos(0) l+cos 2 (fl) 



r + iA 13 - ^(cos(2t9)A 13 - sm(26)p y ) 
ii+co B (29) ro + _ COS (20))A + sm(20) Py 



(50) 



where the matrices and p y are defined as in Tab. I, and, 
for the sake of orthogonality, A := §1 + A i2 + 5A13. 

Therefore, the Lie wedge can be expanded within the five- 
dimensional Hilbert space Ti. := span{H x , H z ,p y , A, Tq} and 
the final inner approximation to the Lie wedge takes the form 



W 



(51) 



<H Z > 




<H X > 



Figure 2. (Colour online) The Lie wedge of the more general system in 
Example 3 is five-dimensional. Here the surface of the convex cone co 
generated by e^^ro + H z )e~ eH y (see (52) with A = 1) is projected 
into the subspace span{H x , H z , To}. Since [To,H y ] 7^ it deviates from 
the rotational symmetry of Example 2 given in the inner part for comparison, 
cp Fig. 1(a). The figure scales with A — > 00. The prism-shaped projection into 
the subspace sprni{H y , H z , To} (not shown) is similar to that of Example 2 
already given in Fig. 1(b). 



where Co is parameterised (again in the short-hand of (41)) as 



Co := 



2sin(0) 
2cos(0) 
7sin(20) 
7(l-cos(20)) 
(ll+cos(20))/6 





" -Hx" 








Py 




A 




. To . 



e g 



}■ (52) 



It is shown in Fig. 2. — As in Example 2, letting Adcxp^n,) 
act on the drift terms adds no further elements to the edge of 
the wedge, so one gets: 

Proposition 3.3: The set tr>o = M.H y (—Co) is the global 
Lie wedge to control system (45). 

Proof: The Lie wedge property of roo can be derived as 
in Example 1. Then the globality of too follows again from 
Theorem 3. ■ 

Generalising the relaxation operator in Example 3 to To := 
7 diag(a,6, c) with a, 6,0,7 > results in a generalised p- 
component replacing (50) by 



Ry{6)T O Ry(0) 
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a+(c-a) sin 2 (8) (c-a) sin(0) cos(e) 
6 

(c-a) sin(0)cos(0) c-(c-o) sin 2 (6) 



= T + j(c - a) sin 2 (6»)A 13 + j{c - a) sin(0) cos(0)Pj, 



= r n 



-r(c-a) a 
2 " 
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7(c~°) | 



cos(2(9)Ai 3 - sm(29)p y 



This leads to a cone Co for (51) that keeps the structure of (52) 
in a slightly more general form, where Example 2 is readily 
reproduced by c = a, while Example 3 follows for a = b = 1 
and c = 2. 

The Lie wedge in Example 3 and its generalised form 
treated above do not take the form of a Lie semialgebra either. 
Choose B := H y from the wedge too of (51) and recall 
A := r + H z . Then the BCH product 



A*B 



H z — o(H x 



- Py 



(53) 
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leads outside the Lie wedge of (51) since the component H x + 
p y is not within 4 the cone (52). 

As pointed out already, in this section, we have chosen a 
representation in R 3 in order to visualise the Hamiltonian parts 
of the respective quantum dynamics by S*0(3)-rotations. In 
quantum mechanics, this picture can be recovered in the so- 
called coherence-vector representation [27]. Therefore, when 
taking an explicit spin-i representation of ad slt (2) in the 
following chapter, the key results obtained here in Examples 1 
through 3 will show up again. 

IV. Open Single-Qubit Quantum Systems 

In this section, we analyze the standard single-qubit unital 
quantum systems beyond their purely dissipative evolution by 
allowing for Hamiltonian drifts and controls. In view of steer- 
ing open quantum systems, this is an important generalisation. 

A. Markovian Master Equation in Qubit Representation 
Based on the Pauli matrices 



"0 


f 




-i 




'1 " 


1 





, <T y := 


i 


, cr 2 := 


-1 



(54) 



in this section we deliberately depart from the previous nota- 
tion by using the explicit spin-i adjoint representation carrying 
the spin-quantum number j = | as prefactor in given by 

& v ■= |(1 2 ® (?v - crj <8 h) , (55) 

for v £ {x, y, z}, where ® denotes the Kronecker product of 
matrices. One easily recovers the su(2) commutation relations 

[i a p , i a q ] = -e pqr i <5y (56) 

to convince oneself of ad 5U ( 2 ) := {ia x ,ia yi ia z }u e = ad su ( 2 ). 
Here and henceforth we use e pqr to discriminate even and 
odd permutations of (x,y,z) by their signs, i.e. e pqr = +1 if 
(p,q,r) is an even permutation of (x 7 y,z), while e pqr = — 1 
for an odd permutation. 

Thus for a single open qubit system in the above represen- 
tation the controlled master equation (13) or rather its group 
lift (16) takes the explicit form 

X{t) = -(^(Ha + ^UjHj) +f x ) X(t). (57) 

3 

Here, X(t) may be a density operator regarded (via the 
so-called vec -representation 5 or a qubit quantum channel 
represented in GL(A, C). Moreover, Hd and Hj are in general 
of the from Hd := (1 2 <8> Hd — Hj ® 1 2 ) and similar Hj. To 
ensure complete positivity, the relaxation term Tl shall be 
again of Lindblad-Kossakowski form which for the standard 
unital single-qubit systems (with Vk Hermitian) simply reads 
Tl = 2^fc7fc cri to give the nicely structured generator 

C u = i(H d + UjHj) +2 ]T ik o\ . (58) 

j k£{x,y,z] 

4 This would require the equality 2 sin(0) = sin(20) to hold for non- 
trivial 8 beyond its actual solutions 6 = mod n. 

5 Note that the vcc -representation as well as the vector of coherence 
notation just provide different 'coordinates' for the abstract master equation 
(13). [31]) as an element in C 4 



The generator is of this form because the iH terms are in 
the t-part of the Cartan decomposition of g[(4,C) into skew- 
Hermitian (?) and Hermitian (p) matrices, whereas the a\ 
terms are in the p-part. 

B. Single-Qubit Systems Satisfying Condition (H) 

Here we consider the class of fully Hamiltonian controllable 
unital single-qubit systems whose dissipation is governed by 
a single Lindblad operator a\ for some k & {x, y, z} i.e. two 
of the three prefactors jx,Jy,Jz have to vanish. 

Similar to Example 1 of Sec. Ill, choose the controls a x 
and a y to see that such a system fulfills condition (H), since 
{ia x , ia y )u e = ad su ( 2 )- Then it is actually immaterial which 
single Pauli matrix is chosen as the Lindblad operator a\, 
because all of the Pauli matrices are unitarily equivalent. So 
without loss of generality, one may choose k = z, i.e. -y x = 0, 
j y = 0, and j z =: 7. 

Therefore the fully Hamiltonian controllable version of the 
bit-flip, phase-flip, and bit-phase-flip channels are dynamically 
equivalent in as much as they have (up to unitary equivalence) 
a common global Lie wedge 

too := ad 5u(2) © -c , (59) 
where the cone Cn, is defined by 

Co :=R£conv{Ua*Ut \ U £ SU(2)\ (60) 

with 

U:=U®U. (61) 

Clearly, the wedge ron. is global by Corollary 2.1 or, alter- 
natively, by Theorem 3 and its edge -E'(tno) is given by the 
Lie subalgebra ad su ( 2 ). The above Lie wedge is isomorphic 
to the one in Example 1 of Sec. Ill for the particular choice 
that r = diag(l, 1,0). 

C. Single-Qubit Systems Satisfying Condition (WH): One 
Lindblad Operator 

Here we discuss an important class of standard single-qubit 
systems which are particularly simple in three regards 

(i) their dissipative term is governed by a single Lindblad 
operator, Tq := 2"~fd~l for some k G {x,y,z}; 

(ii) their switchable Hamiltonian control is brought about by 
a single Hamiltonian a c for some c 6 {x, y, z}\ 

(iii) their non-switchable Hamiltonian drift is &d for some 
d e {x,y, z}. 

Applying the algorithm for the inner approximation of the 
Lie wedge, we get in step (1) 

tt> c dk (l) := i Ra c © -R+(ia d + 2j&l) , (62) 

where again we note the separation by t-p components. In step 
(2) we identify the span generated by the control ia c as the 
edge E(vo) of the wedge. So the conjugation has to be by the 
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-i29a c 



o+ivcr' 



® e" 



Thus in 



control subgroup, i.e. by e 
step (3) one obtains as t-component of the conjugated drift 

K c d {6)-= e- ie& e(ia d )e ie& ° 

{i &d for c = d (63) 

i cos(d)crd + i £cdq sm(6)fr q else 
and as p -component 



Pk(0) 



= (2 7 a 



2\ i8a c 



2j (cos (9)fr k + £ ckr sm(6)a r y 



for c 
else . 



(64) 



The last expression (for c ^ fc) can be further resolved using 
the anticommutator {A, B} + := AB + BA 



" cos 2 (8) ' 




sin 2 (6) 




. cos(6) sin(e) . 





£cfcr{<5'fc,Oy} + 



2 

l+cos(20) 




1 


l-cos(20) 






sin(20) 




_ -ecfcr((o-J®(Tr) + (c r J(8lO-fe)) _ 



(65) 

where the latter identity gives a decomposition into mutually 
orthogonal Pauli-basis elements. 

To summarize, if the control Hamiltonian neither commutes 
with the Hamiltonian part nor with the dissipative part of the 
drift, one obtains in terms of the above K^{9) and P k {9) 



>-dk ■- 



+ conv{jq(0)+P fc c (0) 



} 



(66) 



However, if [cx c ,r ] 
(66) simplifies by P fe c (6>) 

Cdfc =R(f convj 



0, then the convex cone in equation 

r to 



"cos(0)~ 




i&d 






sin(0) 








9 e 


1 




To 







(67) 



in entire analogy to Example 2 of Sec. III. 

The final Lie wedge admits the orthogonal decomposition 



dk 



'■ilk 



(68) 



and moreover by Theorem 3 (or alternatively by Corollary 2.2) 
it is global. For 7 > 0, the edge E(to) = (i(J c ) is again the 
span generated by the control, yet it flips into the full algebra 
E(m) = ad su ( 2 ) m me limit 7 = 0. 

The relation to Examples 2 and 3 of Sec. Ill is obvious: 
Let a unital qubit system satisfy the (WH)-condition and have 
a dissipative Lindbladian To := 270^, induced by a single 
Lindblad operator \% = <Jk- If [coTo] 7^ 0, one arrives at a 
situation resembling Example 3, whereas if [<r c , T ] = 0, one 
obtains a result analogous to Example 2. 

Application: Bit-Flip and Phase-Flip Channels 

Also the relation to standard unital qubit channels is imme- 
diate: Note that in the bit-flip channel the noise is generated 
by <7^, while it is fry in the bit-phase-flip channel and fr 2 z in the 
phase-flip channel, see Tab. II. In the absence of any coherent 
drift or control brought about by the respective Hamiltonians 



Hd = frd or Hj = frj, the Kraus representations are standard. 
By allowing for drifts and controls, the Kraus rank K of 
the channel usually increases to K=A with exception of a 
single frd or frj commuting with the single Lindblad operator 
fr\ keeping K=2. Also the time dependences become more 
involved. Hence explicit results will be given elsewhere. 

Under full H-controllability, the Lie wedges of all the three 
channels become equivalent as the Pauli matrices and thus the 
corresponding noise generators are unitarily similar. 

In contrast, for the case satisfying the (WH)-condition, 
assume a control system with a Hamiltonian drift term gov- 
erned by fr z . Upon including relaxation, now there are two 
different scenarios: if the control Hamiltonian (indexed by 
c G {x,y, z}) commutes with the noise generator (indexed 
by k £ {x,y, z}), one finds a situation as in Example 2 and 
(67), otherwise the scenario is more general as in (66). 

D. Single-Qubit Systems Satisfying Condition (WH): Several 
Lindblad Operators 

Consider a unital qubit system satisfying the (WH)- 
condition and whose Lindbladian To is generated by I = 2 
or I = 3 different Lindblad operators fr\. Then one obtains 
the following generalisations of the symmetric component 



For 1 = 2 and a c _L <jfc, a, 
PkA0) = ' 



1 

7cos 2 (0) 
7 sin 2 (6>) 
7 cos(0) sin(0) 



7 

2(7+7') 
7 (l+cos(2e)) 
7(l-cos(2S)) 

7sin(20) 



\ 

-(o-J(g)o- fc ) 
-(ovT®oy) 

cfcr((o-J«KTfe) + (<Tfc 



(69) 



while for I = 3 and a c _L cr^., a c _L ov, <r c 



kk' k 1 



•(f) 



7 

7COS 2 (0)+7' sin 2 (0) 
7' cos 2 (e)+7sin 2 (0) 
(7-7') cos(6»)sin(e) 



Ok" 



2(7+7'+7") 
7+7'+(7~ 7') cos(28) 
7+7' — (7— 7') cos(2S) 
(7-7') sin(20) 



1 

-(Ofc" ®CT fe ) 
-(<r£®<7 fc ,) 



which for 7 = 7' = 7" simplifies to 

P c kk , k „{6)=T = 2 1 {frl + frl,+frl„) 



(70) 



(71) 



Note that (69) with 7 = 7' precisely corresponds to Example 3 
in Sec. III. 



Application: Depolarising Channel 

Treating the depolarising channel also becomes immediate, 
since one has three noise generators governed by all of fr x , fr y , 
and fr z . Thus the fully Hamiltonian controllable version of 
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Table II 

Controlled Single-Qubit Channels and Their Lie Wedges 



Channel 


Primary* Lindblad Operators 


Primary* Kraus Operators 


Lie Wedges 














case satisfying (WH)-condition 


H-controllable case 


Bit Flip 


Vi = , 


/(111 <?x 


Ei = 




»L = (**«> © 


tt> = ad su(2) e -c 








E = 




[see Eqns. (66,67)] 


[see Eqn. (60)] 


Phase Flip 


Vi = , 


/022 cr z 


Ei = 




X0 dz - (lcr c ) -C dz 


— same as above — 








E = 




[see Eqns. (66,67)] 




DIL-rndSc Flip 


V 1 = . 


/O33 


E x = 


y/^33 <?y 




— same as above — 








E = 




[see Eqns. (66,67)] 




Depolarizing 


Vi = , 


/(ill fa 


Ei = 




d,xyz \ L / ^ d,xyz 


IT) = ad BU ( 2 ) ffi -Crri/z 




V2 = 


/Q22 cr y 


E 2 = 




[see Eqn. (63) and Eqns. (70,71)] 


[see Eqns. (72,73)] 




Vs = 


/033 °"z 


E3 = 


y/T3 cr z 












Eq = 













*) Primary operators are for purely dissipative time evolutions (no Hamiltonian drift no control). Then the time dependence of the 
Kraus operators roots in the GKS matrix {aa}^^. Define: Ai := a22 + Q33, A2 '■= 0,22 — 133, A3 := an + 022, and thereby qu := ^ (1 + e~ aii± ), 
ni ■= 1(1 - e - a «*) and r := \(1 + e" x ^ + e~ A 2* + e~ A 3*), ri := \(1 - e - *! 4 + e -A 3* - e~ A 3*), r 2 := |(1 + e- A i' - e~ A 2* - e - A 3*), 

r # 3 := 1(1 — e~ Alt — e~ A2 ' + e - A3 '). — Under Hamiltonian drift and control the Kraus-rank gets K = 4 except for one single control or drift f/j 
that commutes with the only Lindblad operator V%: in this case the Kraus rank is K = 2. Time-dependences are involved and will be given elsewhere. 
The Pauli matrices cr„ are defined in Eqn. (54). 



the depolarising channel follows the bit-flip and phase-flip 
channels in the structure of its global Lie wedge 

ro := adsii( 2 ) © ~^x V z , (72) 
where the cone c xyz now reads 

c xyz := R+conv{U( lx a 2 x + ly a 2 y + lz a 2 z )U^ | U e SU(2)} 

(73) 

with U of the form (61). Again, the edge of the wedge is given 
by the entire algebra -E'(tn) = ad su ( 2 ) and globality of the 
wedge follows by Theorem 3 or Corollary 2.1. — Moreover, 
note that the Lie wedge in the fully Hamiltonian controllable 
depolarising channel with isotropic noise takes the structure of 
a Lie semialgebra as (in the coherence-vector representation) 
it corresponds to the special case of Example 1 in Sec. Ill, 
where the relaxation operator is a scalar multiple of the unity, 
Tq = A • 1. For anisotropic relaxation, however, this feature 
does not arise. 

If only condition (WH) is satisfied, there are two distinc- 
tions: if the noise contributions are isotropic (i.e. with equal 
contribution by all the Paulis through j x = j y = j z ), one finds 
a cone expressed by (63) and (71). However, in the generic 
anisotropic case, the cone can be expressed by (63) and (70), 
see also Tab. II. 

V. Open Two-Qubit Quantum Systems 

In this section we extend the notions introduced in the 
previous chapter to three types of two-qubit quantum systems 
with uncorrected noise. The two qubits will be denoted 
A and B, respectively. Moreover, we use the short-hands 



'■= S5 0V with p, v E {x, y, z, 1}, where o\ := 1 
as well as the corresponding 'commutator superoperators' 

o> := |(1 <g> o> - r7*„ (g> 1). 

A. Fully H-Controllable Two-Qubit Channels 

A fully Hamiltonian controllable two-qubit toy-model sys- 
tem with switchable Ising-coupling is given by the master 
equation 

P = ~ {i ^ UjZj + Tp) p (74) 
j 

where &j G {&xu o y \\ &i x , o\ y \ a zz } are the Hamiltonian 
control terms with amplitudes {xij} J=1 6 M. 

Since (iaj \j = 1,2,..., 5) Lie = ad su ( 4 ), the edge of the 
wedge is E(n>) = ad 5u (4). Following the algorithm for an 
inner approximation of the Lie wedge, step (1) thus gives 

Wi := ad su(4) © -K+r . (75) 

Conjugating the dissipative component by the exponential map 
of the edge and then taking the convex hull yields the convex 
cone 

c ^convjAcTroC^lC/ :=U<Z)U,U G SU(A),X > 0} , (76) 

which is the two-qubit analogue of the cone in Eqn. (60). The 
resulting Lie wedge 

Itjo := ad su(4 ) © -co (77) 
is global by Theorem 3. 
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B. Two-Qubit Channels Satisfying the (H)-Condition Locally 
and the (WH)-Condition Globally 

By shifting the Ising coupling term from the set of switch- 
able control Hamiltonians into the (non-switchable) drift term, 
&d = &zz, one obtains the realistic and actually widely 
occuring type of system 



(78) 



where now one just has the local control terms &j 6 

{CTasijCTyijoia^o-ij/}. Since (icr x i,i<J y i)\_\e = ad SUA (2^<g> 1b, 
whereas on the other hand {id\ x , iai y )ue = 1a <8> ad SUB ( 2 ), 
the edge of the wedge 

E(H>) = ad 5UA(2)§SUB(2) (79) 

is in fact brought about by the Kronecker sum of local algebras 

sua (2) ® 1b + 1 A ® su B (2) =: su A (2)®su B (2) (80) 

forming the generator of the group of local unitary actions 

exp (sua(2)§0u b (2)) = SU A (2) ® SU B (2) . (81) 

Remarkably, in this important class of open quantum- 
dynamical systems, qubits A and B are locally 
(H)-controllable, respectively, while globally the system 
satisfies but the (WH)-condition. 



The final Lie wedge in these systems reads as 



<k = ad * 
with the convex cone 



_ r 202 

u A (2)©su B (2) w t dfe 



<T ■= K 



conv + J*® 8 } 



(82) 



(83) 



being given in terms of the respective t and p-components. 
Here we use the short-hand of (61) in the sense of C/ 2 ® 2 := 
U 2 (g,2 ® U 2 ®2 to arrive at 



C. Two-Qubit Channels Satisfying Only the (WH)-Condition 

In the final example of a two-qubit system, the independent 
local controls shall even be limited to either x or y-controls 
on the two qubits according to 



p = -(i(&d + UAO-cl + UBCTlc') + T )p 



(88) 



where now &d :— i(& z i + &\ z + a zz ) and a c \ with a single 
c G {x, y} and likewise &i c > with a single d G {x, y} and 
ua,ub £ R. Furthermore, assume the system undergoes local 
uncorrelated noise in each of the two subsystems in the sense 
that the Lindblad operators are of local form 

V k e {ff fc i | k e {x,y,z}} (89) 
V k , G {* lfc , |fc'G{x,j/,z}}, (90) 

where and k' are chosen independently k, k' G {x, y, z} so 
that in the convention of (55) one finds 



T := + 2ia\ k , 



(91) 



This system satisfies but the (WH)-condition both locally and 
globally, the latter following from 



(icr c i,iai c ',i&d)ue = &4»u(4) . 



The Lie wedge is given by 

tofcfc' := (i°c) + (ted) 



'-kk' i 



(92) 



(93) 



where the two-dimensional edge of the wedge is generated by 
the rays (ia c i), (io'ic 1 ) an d the cone 



C c k i := K+ conv {K c {6) + K c ' {6') + K cc '{9, 9') 



(94) 



is given in terms of the t- and p-components (setting 6 := it a 
and 9' := ub and using the relations in (63)) as 



K c (8)+K c (9')+K cc (6,9') 



cos(e) 






sin (9) 






cos(e') 






sin(e') 


■i 




cos(8) cos(# ) 






cos(0) sin(fl') 






cos(e') sin(S) 






- sin(fl) sin(e') - 







(95) 



AM 



{EW^d)^ | [/ 2 ® 2 G 5t/(2) (8) S*7(2)} (84) and (as in (65)) 

t . 

P C k(0) = 27 



t® 2 : = {tWr o )!7 2 T 02 I U 2m G 5/7(2) ® Stf(2)} . (85) 



As before, this immediately results from the initial wedge 
approximation by step (1) 



" cos 2 (fl) " 




sin 2 (0) 




_cos(0)sin(0)_ 





Eckr{&kl ,6-rl } + 



* 202 ~ ad 5UA(2)&UB(2) 9 -MZ(ifid + To) (86) 
followed by conjugation with Ad exp £( ro ) = Ad 2l g, 2 to give 

(87) 



2 




1 


l+cos(2fl) 




-(o-jJ"(8<7h)8lB 


l-cos(20) 




-(ct-J®oV)®1b 


sin(20) 







(96) 



K 2®2 + p2®2 0su(msu(2) (iu d + To) 



as well as 

P$(8') = 2 1 ' 



cos 2 (8') 
sin 2 (9') 
cos(S') sin(fl') 



r ~2 -2 i 



Step (3) then takes the convex hull. — To show globality, let 
to' denote the global Lie wedge corresponding to the fully 
H-controllable system given in Eqn. (74). Then to 2 ® 2 C to', 
and it can be shown that tn 2 ,® 2 satisfies the conditions of 
Corollary 2.2 and therefore is global. 



2 




1 


l+cos(20') 






1— cos(2fl') 






sin(26>') 




_ -e c , k j r , l A ®((o-J,®o- r /)+(o-J,®o- fe /)) 



(97) 
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To see this, observe that by step (1), the initial wedge 
approximation is given by 

tti := (ia e ) + (i&c) © -M^(ia d + 2^1 + 2ia 2 k ,) , (98) 

which has to be conjugated by Ad exp ^E(m))- As usual, the 
edge of the wedge is invariant under such a conjugation, so 
we need only determine the effects on the drift components of 
the system as is done in Eqns. (95) through (97). Moreover, 
the wedge is global by application of Corollary 2.2. 

Now, the generalisation to systems with more than two 
qubits satisfying the (H)- or (WH)-condition is obvious: 
assuming uncorrelated noise, the p-parts of the Lie wedges 
can be immediately extended on the grounds of the previous 
description, since all processes are local on each qubit. Though 
straightforward, calculating the ^-components becomes a bit 
more tedious: but the many-body coherences have to be 
considered just as in (95). 

VI. Outlook: Approximating Reachable Sets 

Knowing the global Lie wedge of a coherently controlled 
Markovian system provides a convenient means to efficiently 
approximate its reachable sets. As in the case of a Lie algebra, 
the image of the wedge ro under the exponential map yields 
a first approximation of the corresponding Lie semigroup S. 
Unfortunately, this image is in general only a proper subset 
of S — this, however, may happen also for Lie algebras when 
the corresponding Lie group is non-compact. Therefore, one 
has to allow for finite products of the form e Al e A2 ■ ■ -e Al 
with Ai,A2,...,A( G ro to obtain the entire semigroup 
S. Although the minimal number of factors to generate 
S (called number of intrinsic control-switches) is in general 
unknown, this approach provides a much more effective 
parametrization of the reachable sets than the standard method 
which works with the original control directions and piecewise 
constant controls as parameter space. Thereby one can opti- 
mize target functions almost directly over the reachable sets 
thus complementing standard optimal control methods of open 
systems [32-34]. Particularly simple are systems whose Lie 
wedges do carry a Lie-semialgebra structure (like in isotropoic 
depolarising channels). Here one knows a priori that only a 
few (or sometimes even zero) intrinsic control -switches are 
necessary, so some control problems may actually be solved 
by constant controls. 

VII. Conclusions 

We have generalised standard unital quantum channels (bit- 
flip, phase-flip, bit-phase-flip, and depolarising) by allowing 
for different degree of coherent Hamiltonian control. For 
the first time, here we have characterized their respective 
global Lie wedges governing all directions the controlled open 
system can possibly take. The results have been further gener- 
alised to various types of two-qubit systems with uncorrelated 
noise. Since controlled multi-qubit channels can be treated 
likewise, the geometrical Lie-semigroup approach taken is 
anticipated to find wide applications in quantum systems 
theory and engineering: this is because knowing the global 



Lie wedge of a controlled Markovian system paves the way to 
efficiently approximate its reachable sets. Thus this knowledge 
will be very useful for improving known bounds (cf. [16]) on 
the corresponding system semigroup Ps in follow-up work. 

Finally, our results demonstrate that the Lie wedges asso- 
ciated to most of the controlled quantum systems do not take 
the special form of Lie semialgebras, an important exception 
being the fully controlled isotropic depolarising channel. 

VIII. Appendix 
A. The Principal Theorem of Globality 

For the reader's convenience, we state the 'Principal Glob- 
ality Theorem' with minor simplifications. For the full version 
and its (quite involved) proof we refer to [2], which we sketch 
in the sequel. 

Let G be a matrix Lie group with Lie algebra g, so 

S X:={AX\Ae 9 } (99) 

can be envisaged as tangent space TxG at X G G, while TG 
and T*G shall denote the tangent bundle and, respectively, 
cotangent bundle of G. Thus, one has the isomorphisms 

TG = g x G and T*G g* x G. (100) 

Now, let ro be any wedge of g. A 1-form on G is a smooth 
cross section of the cotangent bundle, i.e. u : G — > T*G with 
oj(X) G T£G. Moreover, uj is called 

(1) exact if there exists a smooth function Lp : G — > K such 
that d(p = w; 

(2) to-positive at X <= G if (u(X),AX) > for all A G ro; 

(3) strictly to-positive at X G G if ro-positivity holds at X G 
G and one has (lo(X), AX) > for all A G to \ -ro; 

The existence of a strictly to-positive 1-form is ensured in the 
following scenario [2]: If G is a Lie group with Lie algebra 
g and H a closed subgroup with Lie algebra t), then for any 
Lie wedge ro G g whose edge E(ro) coincides with f) one can 
construct strictly ro-positive 1 -forms on G. Note, however, that 
these 1 -forms on G are in general not exact. Yet, whenever 
exactness can be guaranteed in addition, one has the following 
equivalences. 

Theorem 4 ([2]): Let G denote a finite-dimensional real 
matrix Lie group with Lie algebra g and let ro be a Lie 
wedge of g. Moreover, let go := (to)Lic be the Lie subalgebra 
generated by ro and let Go be the corresponding Lie subgroup 
of G. Further, assume that Go is closed within G. Then the 
following statements are equivalent: 

(a) to is global in G. 

(b) ro is global in Go. 

(c) There is a closed connected subgroup H of Go with 
i(H) = E(tv) and a 1-from u on Go which satisfies 
the following conditions: 

(i) ui is exact. 

(ii) uj is to-positive for all X G Go- 

(hi) iu is strictly ro-positive at the identity 1. 

Now, the following consequence of the 'Principal Globality 
Theorem' already mentioned in the main text is a useful tool 
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whenever a global Lie wedge to embracing the Lie wedge too 
of interest is already known. 

Corollary 2.2 ([2]) Let G be a Lie group with Lie algebra 
and let too C ro be two Lie wedges in a. Provided 



ro \ —too c ro \ — ro , 



(101) 



then too is global in G if the following conditions are satisfied: 

(i) tt> is global in G. 

(ii) The edge of too is the Lie algebra of a closed Lie 
subgroup of G. 

In other words, if the edge of the wedge follows the 
intersection -E'(rDo) = E(to) n too and ro is global, then too 
is also a global Lie wedge, whenever cxpS(too) generates a 
closed subgroup. 

B. Lie Semialgebra Structure in Example 1 

In Sec. Ill we stated that the Lie wedge of Example 1 

to :=so(3) © (-c ), (102) 

where c := R%M(T ) and M(T ) := {S G snm(3) | S -< 
To}, is in fact a Lie semialgebra for To = X- 1 (corresponding 
to the isotropic depolarising channel), whereas it fails to be 
a Lie semialgebra for any other To- Recall, here stjm(3) is 
the set of all symmetric 3x3-matrices. For proving the above 
statement, we distinguish the following cases 6 : 

(i) To is a multiple of the identity, thus we can assume 
without loss of generality To := 1; 

(ii) To has zero as eigenvalue with multiplicity 2, thus we 
can assume To := diag (1, 0, 0). 

(iii) To has an eigenvalue different to zero with multiplicity 
2, thus without loss of generality To := diag(l, 1, 0). 

(iv) r has three distinct eigenvalues, i.e. r :— diag(a, b, c) 
with a > b > c > 0. 

In all cases, the identification of the dual wedge of too 
is crucial for the compution of the tangent space T^tUo at 
A G too via (6). Therefore, we first provide an auxiliary result 
characterizing the dual cone of Co within srjm(3). 

Lemma 8.1: Let To := diag (a, b, c) with a > b > c > 
and let c := R£M(ro) with M(T ) := {S G srjm(3) | S < 
To}. Then the dual cone of Co within srjm(3) is given by 

cj := {^GSt 5 m(3)|cA 1 (5) + feA 2 (5)+aA 3 (5) >0},(103) 

provided Ai(S) > A2(5) > A3 (S) are the eigenvalues of S. 

Proof: By definition, one has the equivalence S G 
c* if and only if (S,S') > for all S' G c . 
Since Co = R^ conv Oso(3) C^o) this condition reduces to 
(S,er e T ) > for all e e SO(3). Then von Neumann's 
inequality [35] provides the equivalence: (S, QT Q T ) > for 
all 9 G SO(3) if and only if cXx(S) + b\ 2 (S) + a\ 3 (S) > 0, 
where Xi(S) > X2(S) > X 3 (S) are the eigenvalues of S. 
Hence the result follows. ■ 
Now, we are prepared to prove the above claim about the 
Lie semialgebra property of too 

6 Although case (iii) seems to be quite similar to case (ii), its proof is more 
involved and a helpful preparation of the general case (iv). 



Proof: (i) In case T = 1, the pointed cone Co := 
RQ.M(r ) equals the ray Rp 1. By Lemma 8.1, we obtain 
c; = {S G st)m(3) I trS> 0} and thus = sl(3, M) ® K+r , 
where sl(3,R) denotes the set of all 3x3-matrices with trace 
zero. Hence 

too = so(3)- L n Cq = {S G snm(3) | tr S > 0} . (104) 
For A := Al + O G to with A > and VL G So(3) it follows 
{S G snm(3) | tr S = 0} for A > 0, 



A 1 - n too = 



and thus 

T A to = (A 1 - n to*) 



{S G snm(3) | tr S > 0} for A = 0, 
so(3)©Rl for A > 0, 



50(3) 



for A = 0. 



(105) 



(106) 



Thereby the inclusion [A, T^too] C T^too is obviously always 
satisfied and hence too is a Lie semialgebra for To = 1. 

(ii) In case To = diag (1,0,0), it is easy to see that the 
pointed cone Co := R qM(To) actually consists of all positive 
semidefinite 3 x 3-matrices. Here, Lemma 8.1 yields c*, = Co. 
This reflects the well-known fact that the cone of all positive 
semidefinite matrices is self-dual within the space of all 
symmetric matrices. Hence 



to. 



Now, for A 



so(3) x n c* = snm(3) n Co = c . (107) 

G too we obtain 



"100" 




'0 


-1 


0" 





+ 


1 



















0. 



A n too = 

and therefore 
T A to = (A 1 - n ro* 





s 



S G snm(2) , S > 



(108) 



00(3) ® span{T , p y ,p z }. (109) 



Finally, for disproving the inclusion [A, T^too] C T^too 
consider the commutator of A G too and B := p z G TaWo • It 
follows [A, B] = — H z + diag (—2, 2, 0) which clearly violates 
the inclusion [A, T^too] C T^too- Thus too is not a Lie 
semialgebra for To = diag (1,0,0). 

(iii) In case To = diag (1, 1,0), we obtain by Lemma 8.1 the 
following description 

ro* Q = so(3) ± nc^ = cj = {S e snm(3) | A 2 (5)+A 3 (5) > 0} . 
Now, let A := To + H y . Then, it is easy to see that 
A x nto* D 

R q conv {diag (1, -1, 1), diag (-1, 1, 1), (p z + E 33 )} . 
Moreover, for S G A n top one has the conditions 

(r ,5) = o and (er e T ,5)>o 

for all Q G SO(3). Now, differentiating the second condition 
with respect to 9 G SO (3) shows ([fi,r ],S) = for all 
G so(3), i.e. [fi, T ] belongs to {A 1 - n too)" 1 ". Thus one has 

T^too = (A^ntoo)^ 3so(3)©RA©span{p :c ,py}. 
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Hence, counting dimensions finally yields 

T A w =so(3)®span{T ,p x ,p y } . 

To disprove the set inclusion [A, T^too] C T^too consider 
the commutator of A G too and B := p y G TaK>o . The 
computation is left to the reader (see Tab. I). The result clearly 
violates the inclusion [A, Ta too] C T A K>o an d thus too is not 
a Lie semialgebm for Tq = diag (1, 1,0) either. 

(iv) For Tq = diag (a, b, c) with a > b > c > and A := 
To + H v with v G {x, y, z}, the same arguments as above 
show that TaK>o is given by 

T A tv = so(3) ®span{T ,p x ,Py 7 p z } ■ 

Therefore, an appropriate choice of B = p v with v G {x, y, z} 
demostrates again that too is not a Lie semialgebm in the 
general case To = diag (a, b, c) with a > b > c > either. ■ 
Note that in all the above cases the tangent space of too has 
the following form 

TAtoo =so(3)®Rr o ©Tr o 0so(3)(ro), 

where the tangent space of the orbit Oso(3) (ro) at Tq is given 
by Tr o S o(3)(ro) = m, To]\ne so(3)}. 
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